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Abstract In the presented study, a numerical model which predicts the flow-induced collapse within the 
pharyngeal airway is validated using in-vitro measurements. Theoretical simplifications were considered to 
limit the computation time. Systematic comparisons between simulations and measurements were performed 
on an in-vitro replica, which reflects asymmetries of the geometry and of the tissue properties at the base of 
the tongue and in pathological conditions (strong initial obstruction). First, partial obstruction is observed and 
predicted. Moreover, the prediction accuracy of the numerical model is of 4.2 % concerning the deformation 
(mean quadratic error on the constriction area) . It shows the ability of the assumptions and method to predict 
accurately and quickly a fluid-structure interaction. 

Keywords numerical simulation and modelling • in-vitro measurements • fluid-structure interaction ■ 
obstructive sleep apnea syndrome 

1 Introduction 

Since the 90's, biomechanical modelling of the human upper airway has received a growing interest since it 
allows a better understanding of its physiology and pathophysiology, as well as an increased quality of the 
treatments of its specific pathologies. Among these, Obstructive Sleep Apnea Syndrome (OSAS) has been 
object of particular attention, as it became a major health care topic, affecting a growing part of the pop- 
ulation, especifically in Europe and in the United States [22jf39] . It is characterized by the occurence of an 
abnormal rate of apneas and hypopneas during sleep 22J. During an episode of obstructive apnea (respec- 
tively hypopnea), the soft tissue in the pharynx completely (respectively partially) collapses in response to 
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inspiratory airflow. It induces a temporary cessation (respectively limitation) of the respiration [3j. The main 
health effects are excessive daytime sleepiness and an increased risk of cardiovascular diseases [22] • In parallel 
to the great amount of medical/biomedical engineering research that has been carried out to understand this 
highly complex phenomenon (see e.g. 38, 1 for recent works or 22,3, 12,9, 10 for recent overviews), theoretical 
biomechanical models have been proposed such as simplified mathematical models of the interaction between 
the airflow and the soft tissue [TBTfTlE ]. numerical models of the respiratory fluid flow 28,19 and numerical 
models of the fluid-structure interaction [231121(137] . They allow indeed to understand the relationship between 
the fluid-structure interaction in the pharynx and the abnormal flow patterns observed in the medical or phys- 
iological experiments. They might allow furthermore to predict these flow patterns using measurements of the 
biomechanical properties of the upper airway (geometry, rheology). This makes them of interest to improve 
the quality of the treatments such as surgical procedures [8|fT4] or mandibular advancement splints [15] . 

In a previous paper [6], a numerical model which aims at predicting the fluid-structure interaction within 
the pharyngeal airway has been described. It is based on simplified assumptions so that the computational 
cost of simulations be compatible with clinical applications. An original in-vitro tongue replica has also been 
described. It takes into account the very specific properties of the upper airway. In particular, the asymmetry 
of the pharyngeal duct at this level has been reproduced in a simplified manner. A few comparisons between 
predictions of the numerical model and measurements of the flow-induced collapse within the in-vitro replica 
provided encouraging results (see [3] for details). Nevertheless, this was not sufficent to be really conclusive 
about the interest and drawbacks of the numerical model. Therefore, the aim of this study was extensive 
validation of this numerical model through in-vitro experiments, and to our knowledge, it has not been done 
so far in this context (for instance in [37], the validation is carried out with rigid walls). The main novelty 
in comparison to the study presented in [6] is that systematic comparisons have been carried out, for a wide 
range of parameter values. Moreover, pressure sensors have been integrated into the experimental setup. As 
the pressure at the level of the constriction plays a critical role in the flow- induced collapse (see [34] ), it was 
jugded interesting to assess its prediction in conditions of fluid-structure interaction. Besides, a digital camera 
quantifies the two-dimensional deformation of the simplified airway at the main site of obstruction. This gives 
more detailed and pertinent information than the laser beam which was of use in [6]. 



2 Material and methods 

2.1 In-vitro setup 

The aim of the experimental setup described here is to reproduce a fluid-structure interaction in conditions of 
strong obstruction and quasi-steady motion. It obviously simplifies the complex in-vivo reality, but allows to 
obtain reliable dynamical measurements of the pressure and of the deformation in controlled conditions, which 
is required for validation. 

The in-vitro setup is depicted Figures [TJ & [2j The real morphology of the upper airway has been simplified as 
suggested Figure [TJ (a). However, its overall configuration and the main dimensions have been conserved. In 
particular, a frequent site of collapse is the base of the tongue [25]. This anatomical entity is represented in the 
setup by a short latex tube filled with water (deformable tongue replica). This deformable structure intersects 
orthogonally a rigid metallic pipe which stands for the pharynx. The geometry of the duct at the level of this 
junction is those of a constriction, in accordance with observations at the tongue base. 

The diameter D of the deformable tongue replica is 49 mm. The latex wall (Piercan Ltd.) has a thickness e of 
0.3 mm. Its Poisson's ratio v is of 0.5 (incompressible material). The value of its Young modulus E has been 
determined to 1.68 MPa (see [I] for details). The latex sheet is mounted onto a rigid metallic support (Figure 
[2] (a)). A hole in the support enables the deformation of the latex in response to the fluid flow: the deformable 
part of the cylinder is depicted Figures [2] (a) &Q](b) in dark shade. An external water supply (Figure [21 (c)) 
through a duct connected to a water column allows to control and to measure the water pressure P ext inside 
the latex interface (external flow). Indeed, the height of the water column is controllable [55] . Manipulation of 
this pressure influences simultaneously two parameters: the geometry of the duct and the overall stiffness of 
the deformable part. 
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The rigid metallic pipe in which the airflow (internal flow) circulates is depicted Figure [2] (b). Its diameter d is 
of 25 mm. A removable flat plate allows to change the geometrical characteristics of the duct (Figure [1] (b)). 
Different heights^ of the duct, corresponding to different pharyngeal calibers, can be chosen using different 
plates. Air supply comes from a pressure reservoir, which can be considered as an artificial lung. This is a 
rectangular box of approximatively 0.75 m 3 , fed by a compressor. The control of the pressure in the reservoir 
is ensured by a pressure regulator (Norgren(TM) type 11-818-987). A constant selected pressure, within the 
range 0-3000 Pa, is thus obtained at the inlet of the duct . The pressure within the airflow is measured with 
sensors located at different positions (Figure [T] (c)) One sensor (po) allows the measurement of the upstream 
pressure. Two others {pi,P2) are located at the bottom of the rigid plate, and allow measurements of the 
pressure at the level of the constriction and downstream. The pressure sensors are within holes of diameter 
0.4 mm. They are piezo-electric0 and deliver a tension proportional to the pressure. Preliminary calibration, 
with a liquid manometer (Kimo(TM)), allows to determine the constant of proportionality of the relationship. 
Then, the pressure is measured with a precision of ± 1 mmH^O. 

An important geometrical parameter is the minimal height of the duct, at the bottom of the upper cylinder, 
called the constriction height h c (Figures Q] (b) (c) & 121(c)). As explained before, this height could be changed 
using different metallic plates, while P ex t is maintained. In physiological terms, it corresponds to the minimal 
pharyngeal caliber. Another important parameter is the area at the level of the constriction A c . It corresponds 
to the minimal pharyngeal area. We chose this parameter for the measurement of the flow-induced deformation 
since it captures its global effect on the obstruction of the duct (h c is only a punctual measurement). To this 
purpose, a digital camera has been used (Figure [T] (c)), which measures the bidimensional geometry at the 
level of the constriction, in frontal view and at the downstream end (Figures [T] (b) &[H (c)). The resulting 
information is the constriction height h c (z) for all z, from which we obtain the constriction area A c and follow 
its variation. 

The digital camera is an industrial modeH, with a resolution of 1280 x 1024 pixels. Its focal has been adjusted 
so as to visualize with the best accuracy the constriction site (focal plane located at the site of the constriction, 
and as wide as the rigid pipe). The zoom objective is such that the number of pixels corresponding to the 
constriction area is maximized. The contrast is enhanced with an ordinary light source upstream the pipe. A 
software^ allows to control the parameters of the camera, to visualize and to store the resulting acquisitions. 
Those can be obtained automatically, at periodic time intervals, in synchronization with the pressure mea- 
surements, thanks to LabView^. A preliminary calibration is necessary, to determine the relationship between 
distances measured on the picture and real distances. The method consists in localizing three points on the 
picture. The distance between these points is known. After the calibration, the horizontal and vertical resolu- 
tions are determined: = r v ~ 6.3.10 -3 cm / pixel. 



2.2 Mechanical theory for fluid-structure interaction 

Numerical computation of a complete fluid-structure interaction problem, i.e. solving the incompressible Navier- 
Stokes equations dynamically coupled with a deformable structure in large deformations, is still challenging 
nowadays, both in terms of stability and convergence towards the solution and in terms of computation cost 
[2^1fT51l3Tl[33] . Therefore, the complexity of the description of the ongoing phenomena is reduced, and simpli- 
ficatory assumptions are stated. The further outlined assumptions are made in agreement with the anatomy 
and the physiology of the upper airway, and they allow to reduce the computation time to a great extent. A 
dctailled description of the biomechanical model and of the numerical method can be found in [6] and [4] . In 
this section, only the main characteristics are given. 

Concerning the structure, the assumptions of a quasi-steady motion, with small deformations and displacements 
have been chosen. The constitutive behavior is given by the Hooke law (linear elasticity). As a result, the 
following equations govern the structure deformation: 

1 dimension in y direction 

2 Endevco 8507C, Kulite XCS-093 

3 Inca 311, Philips 

4 Clicks(TM), version 1.4-0, Philips, Industrial Vision 

5 Lab View 7, National Instruments 
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'V-CT + f = 0, 



< a - (i-2v)(i+v) £ i I + TTU £ > (!) 
i£ = i(Vu + (Vu)*), 

where cr is the Cauchy stress tensor, and f is the vector of external forces. E is the Young's modulus, and 
v is the Poisson's ratio, ej is the trace (first invariant) of the small deformation tensor e. u is the vector 
of displacement. The boundary conditions consist in immobility constraints (displacement u — 0) at sites of 
attachement to rigid parts and imposed external forces in the surface of contact with the airflow. 

The fluid flow is considered as incompressible (Mach number of O(10~ 2 )), laminar (Reynolds number of O(10 3 )), 
stationnary (Strouhal number of O(10 -3 )) and bidimensional [Mj. Since the Reynolds number is of the order 
of 10 3 , we use an asymptotic simplification of the incompressible Navier-Stokes equations: 

' ud s u + vdyU = —d x p + d^u, 

-dyP = 0, (2) 

d x u + d v v = 0, 

(u, v) are the longitudinal (x) and transverse (y) components of the fluid velocity, and p is the pressure |17j . All 
the variables are nondimensional: x = x(hoRe)^ 1 , y — yh$ , p = P(pUq u = uIIq 1 , v = vReU^ 1 ] where ho 
is the transversal dimension of the pharyngeal duct, t/ is the mean longitudinal speed, and Re is the Reynolds 
number (Re = Uaho/is, with v the kinematic viscosity of the air). These equations, called Reduced Navier- 
Stokes / Prandtl (RNSP), allow to take into account the boundary layer formation as well as the separation of 
the fluid after the narrowing of the pharyngeal duct, at the base of the tongue [17]. The boundary conditions 
are the following: no slip on the upper and lower walls and a pressure difference AP imposed between the inlet 
and the outlet. 



2.3 Numerical solving of the equations 

As the problem is considered as quasi-steady, the fluid-structure interaction is solved using a segregative 
method: the equations that govern the fluid and the solid are solved alternatively, inside a global loop. The 
pressure gradient in the fluid is imposed gradually, in st steps : AP\ = 0, . . . , AP st — AP max . For each pressure 
gradient APi at step «, the fluid forces on the wall are first computed and imposed. Then, the wall is deformed, 
which changes the fluid flow domain. The fluid forces need then to be computed again. This is the reason 
why a finite number of iterations it should be imposed at each step, until equilibrium of the wall is reached. 
A convergence criterion cv, which corresponds to the maximal displacement between two iterations, allows to 
ensure that at the end of the it iteration, the displacement of the structure is no more significant. Typically, a 
choice of st = 8 and it = 6 has been found sufficient to ensure the convergence of the algorithm (cv ~ 10 -3 mm). 

The continuum equations of the wall are solved using the finite element method 40 . In the context of small 
deformations and of linear elasticity, the relationship between the nodal displacements and the nodal forces 
is linear. As a result, the precomputation of the inverse of the stiffness matrix [K] is done before the fluid- 
structure interaction loop [7]. It saves computation time as nodal displacements are obtained through simple 
matrix multiplication at each iteration. This preliminary step is achieved using a commercial finite element 
solver (Ansys(TM) Software). The finite element model of the tongue replica, in other terms of the latex tube 
filled with water, is depicted Figure O The mesh is constituted of 150 linear elements of 8 nodes, regularly 
dispatched. The nodes in contact with the hollow metal support (Figured (a)) are immobilized (Figure[3]). The 
boundary conditions are such that a bidimensional model would have not been satisfying. This is the reason 
for the choice of a tridimensional model. The incompressibility of the latex is approximated through the choice 
of a Poisson's ratio v of 0.499. The Young modulus has been fixed to E = 1.68 MPa. The water inside the 
latex tube has been taken into account in the model through the application of constant pressure forces on the 
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surfaces supposed to be in contact with the water. 

The RNSP equations that govern the fluid flow are solved using a finite difference method |16j . which is easy 
to implement and allows fast numerical solving. The grid that is used is of size 2000 x 1000. This is sufficently 
high to capture boundary layer formation and separation, as well as jet formation, and increasing the grid 
resolution above this level does not significantly change the values of the pressure distribution (see [J) and [5] 
for details). As the airflow is modelled using bidimensional equations, the finite element model is divided into 
5 thin slices on which are carried out computations of the fluid forces. A finite element model of 10 slices has 
also been built so as to check the impact of this numerical parameter and no significant difference between 
the two models has been found. The resulting pressure distribution is imposed at the surface of the structure 
in agreement with the principle of virtual works, which provides an accurate approximation of the pressure 
load (see [5] for detailed formulas). The codes for computation of the fluid flow and fluid-structure interaction 
have been written by the authors, using Matlab(TM) interfaced with C language. As a result, the duration of 
the computations for the overall algorithm is typically of the order of 20 minutefl This can be considered as 
reasonable for clinical applications. 

2.4 Methodology for comparisons 

A batch of measurements has been carried out, for which the latex thickness (e = 0.3 mm) and the initial 
constriction height (h® = 1.84 mmj3 remained fix. The following parameters were systematically varied: 

— the external pressure P ex t, 

— the maximal inlet pressure P™ ax . 

The external pressure ranges from 100 to 700 Pa. For each value of P ex t, four measurements have been carried 
out, corresponding to different values of P™ ax : 200, 400, 600 and 800 Pa. Measurements are performed 
according to the following procedure: the inlet pressure P e is gradually increased from Pa up to the required 
value P™ ax . At fixed time steps (0.25 s), a synchronized sampling of the following outputs is done automatically: 

— the upstream pressure P e , the pressure P° at x — mm, and the pressure P c 16 at x — 16 mm, as measured 
respectively by the sensors po,Pi,P2 described in Section [2Tl Note that when P e = Pa, P° corresponds 
to the pressure at the level of the constriction. Yet, during the deformation, the new position of the wall 
can be such that the pressure P® does no longer correspond to the pressure at the constriction. 

— the height of the duct at the level of the constriction : h c (z), using the digital camera described in Section 
12.11 The constriction area A c is determined automatically. 

The experiment is stopped when the value observed for the output does not change anymore. The stationnary 
state is then reached, after approximatively 10 seconds. 

For each value of P ex t, the mean experimental curves P®' exp {P e ) and A exp (P e ) that result from the four 
experiments at different values of P™ ax are computed, with the associated standard deviation. A numerical 
simulation is then carried out. The inlet pressure in the simulation is equal to P e and the outlet pressure is equal 
to P c 16 . The theoretical curves P®' slrn (P e ) and A s ^ m {P e ) are compared to the experimental ones: P^ exp {P e ) 
and A^. xp (P e ). For quantitative comparison we compute first the mean quadratic error, in normalized form, 
which is, for the constriction area: 

(3) 

Here, n is the number of points in the experimental sample, (P*, A^ p {Pl)) the coordinates of the i th experi- 
mental point, (Pg, A s J ,m {PD) the coordinates of the corresponding point on the theoretical curve. The constant 
of normalization A^ is the area at the beginning of the experiment, when no internal airflow is circulating 
(P e = Pa) and after a water pressure P ext has been imposed. For the pressure P®, the quadratic error e(P°) 
is computed identically, except the constant of normalization which is the maximal inlet pressure: p™ ax . A 
maximal error has been computed as well: 

6 Dell Precision 330 (TM) workstation, with Pentium(TM) 4, 2 GHz, 1 GO RAM. 

7 at P ex t = Pa 
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e max (A) = -j^r t max \A^{P\) - Af m (Pl)\ , (4) 

(and an identical definition for e max (P )). Finally, the overall prediction performance is evaluated using the 
coefficient of determination R 2 (0 < R 2 < 1 ): 

i? 2 = l-^, (5) 
a y 

a 2 being the variance of the prediction residuals, and a 2 is the variance of the experimental measurements. 
The closer is R 2 to 1, the best is the prediction of the numerical simulation. As a result, two coefficients of 
determination are computed: R 2 P for the pressure at the level of the (initial) constriction x — and R\ for the 
constriction area. 



3 Results 

The results of the comparisons between predictions and measurements are summarized Table [T] for all the 
investigated values of P ex t, and a typical example of the resulting curves P°(P e ) and A c (P e ) is given Figure [U 
for P ext = 300 Pa. First, it should be noticed that the measurements are strongly repeatable. Table Q] indicates 
a mean value of the standard deviation a exp , computed from the four repeated measurements, of f f .3 Pa for 
the pressure P° and of 0.29 mm 2 for the constriction area A c . Variations of this standard deviation with P e 
can be observed Figure [4] 

First, concerning the pressure P c ° at x — mm, the theoretical and experimental curves are in agreement 
with a mean error of less than 8 %, until P e exceeds 500 Pa. After this value, they split, and the value of the 
computed pressure is always lower than the experimental value (of ~ 30 %, see Figure 0] (a)). Moreover, this 
pressure is both negative and decreasing at the origin, but quickly changes to become positive and increasing. 
This could be explained by the translation of the constriction towards downstream as P e is increased and the 
wall is deformed. As a result, the pressure P° moves away from the value of the pressure at the constriction, 
that remains still negative. The most valuable explaination for the difference observed between the theoretical 
and the experimental values of P® is that the simulation underestimates the displacement of the wall in the x 
direction. Considering the shape of the curve P(x) in the area of the constriction, an error of the order of 1 mm 
in the prediction of the displacement should be followed by an error of the order of 100 Pa in the prediction of 
the pressure. These phenomena lead to a mean quadratic error of 9.8 % and a mean value for R 2 p of 0.64 (see 
Tabled]). 

Then, concerning the constriction area A c , the mean quadratic error is of 4.2 % (R 2 A — 0.71). The theoretical 
and experimental curves are depicted Figured] (b). At the origin (P e = Pa), there is no significant error 
between theory and experiments. This is not surprising since the values of A c for P e = Pa and all the values 
of P ex t have been used to determine a Young modulus that gives the best adequation between simulations and 
data (see [1]). For values of P e between Pa and 200 Pa (0 < Re < 1000), a plateau is observed in the exper- 
imental curves. It is explained by inertial effects in the deformable portion, which arc not taken into account 
in the model because of the assumption of quasi-static deformation. Nevertheless, it is not the major source of 
difference between the theoretical and experimental curves. For 200 < P e < 400 Pa (1000 < Re < 1400), A c is 
a decreasing function of P e . This function, in a first approximation, should be treated as linear. In this range 
of values of P e , the agreement between theory and experiments can be considered as satisfying. The maximal 
amplitude of the closure (AA™ ax / A c (P e = 0)) is of approximatively 18 %, and does not change significantly 
with P ex t- This maximal closure is always obtained for P e ~ 400 Pa. Indeed, in the interval (400 < P e < 700 
Pa, 1400 < Re < 1900), a change is observed: the duct reopens at the level of the constriction. This behaviour 
might be explained by an increase of the downstream pressure P c 16 > ea; p ( U p to ~ 180 Pa, after P e = 400 Pa). 
As p^ 6 ' ex P is used in the simulations to impose the values of the outlet pressure P s , the numerical model is 
able to reproduce this reopening, though the strong non-linearity in the experimental curves is not reproduced 
in the simulations. In the last interval (P e > 700 Pa, Re > 1900), the theoretical and experimental curves are 
diverging. Indeed, after this critical value of Re, self-sustained oscillations in the latex structure are initiated. 
No measurement of A c is possible in this case. As the modelling assumptions include quasi-steadiness, such a 
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behaviour can not be reproducted by the simulations. 



4 Discussion 

First, in the in- vitro experiments, flow- induced obstruction has been observed systematically, with a ratio 
AA™ ax /A c (P e — 0) of approximatively 18 %. This confirms the preliminary experimental results of [5J. This 
effect is due to the internal airflow, which is submitted to acceleration at the level of the constriction and in- 
duces pressure losses (Venturi effect). Since the deformation of the latex wall is governed by the local pressure 
difference between the internal and external flow, these pressure losses result in a decrease of the constriction 
area (partial obstruction of the duct). The values of the couple (P e , P s ) during an experiment are such that in 
the quasi-steady regime (200 < P e < 400 Pa), the expiratory phase is reproduced in a very simplified manner. 
Remembering the analogy between the latex cylinder and the base of the tongue, the airflow which circulates 
in the pipe is similar to an airflow which would go from the hypopharynx to the mouth cavity. P e would then 
be the pressure in the hypopharynx, at the base of the epiglottis, and P s the pressure in the mouth cavity, 
approximatively equal to the atmospheric pressure (as in the simplified models described in [20]l3"6]). Consid- 
ering this analogy as valid, the presented experiments would then reproduce the flow-induced obstruction in 
conditions of expiratory flow, which has been observed in some apneic patients or heavy snorers, and may 
be at the origin of expiratory flow limitation 23, 27, 30, 36 . This analogy, already dicussed in [BJ, is however 
limited, since the experiments do not reproduce with accuracy the complex dynamics of a respiratory cycle 
(see e.g. [38]), and was only focused on the quasi-steady phenomena. The range of values of the upstream 
pressure P e , however, is physiological as it conduces to airflow rates of the order of 10 1/min (see e.g. [32] V 
As a result, an extensive study of the response of the tongue replica to different pressure commands, close to 
some typical physiological or pathophysiological cases, may constitute a first perspective of this work. In par- 
ticular, measurements in conditions closer to inspiration (P e equals to the atmospheric pressure, and P s < P e ) 
would allow to reproduce in a simplified manner what happens during an apneic episode, which is known 
to occur during the inspiratory phase. Of course, the change from expiration to inspiration corresponds to a 
simple change of boundary conditions in the numerical model, and simulations in inspiratory configuration 
have been carried out with success from medical data (sagittal radiographics) |5J. Nevertheless, experiments in 
inspiratory configuration would be of interest, on the one hand, to validate the simulations in this case, and 
on the other, to obtain deformations of the latex structure of larger amplitudes. Furthermore, considering the 
actual geometry of the setup, it is impossible to obtain a complete closure of the channel, which would be of in- 
terest as Obstructive Sleep Apnea is associated to complete closure of the pharynx and complete flow cessation. 

Then, concerning the main point of this study, which is the experimental validation of the numerical model 
described in Section [21 this one can be considered to provide a satisfying first approximation for the prediction 
of the flow-induced collapse measured with the in-vitro setup. Indeed, the main quadratic error for the predic- 
tion of the constriction area is of 4.2 % (R\ = 0.71). This error does not vary significantly with the pressure 
P ex t in the external flow, when P ext < 500 Pa. For P ext — 700 Pa, this error is higher and thus this value 
can be considered as a limit of validity for the numerical model. Moreover, the theoretical and experimental 
curves are in good adequation for the range of values of P e (200-400 Pa) corresponding to quasi-steady be- 
haviour and thus relevant to our study. When the upstream pressure is higher (P e ~ 700 Pa), auto-oscillations 
are observed in the experiments. This behaviour has already been reported on symmetrical |26j as well as on 
asymmetrical geometries [6j. It can be predicted using linear stability theory in association with a simplified 
physical model [26j . In our case and according to the assumption of quasi-steadiness, this phenomenon can not 
be reproduced in the simulations. However, it should be associated to snoring or speech production processes, 
and is therefore not relevant here. The simplified asymptotic theory (RNSP) to model the fluid may be the 
first cause of the differences between numerical predictions and measurements. However, systematic compar- 
isons effectuated with help of a rigid tongue replica proved that this theory provides an accurate prediction 
of the pressure distribution for the geometry we used, even if the assumption of small variations in the axial 
direction necessary to derive the equations is not satisfied [33] . From the earlier study pH] , it resulted also that 
the recirculation effects as well as turbulent effects that may occur after separation of the flow had negligible 
impact on the prediction of the pressure distribution. Moreover, no evidence of turbulence has been found 
in the region upstream the point of flow separation [34] . The assumptions stated for the structure explain 
the other part of the discrepancies between theory and experiments. In particular, the measurement of the 
pressure at the level of the (initial) constriction {x = mm) revealed that the numerical model underesti- 
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mates the displacement of the tongue replica. This might be due to the assumption of small displacements 
which is not respected in this situation. An improvement would then consist in using a shell or membrane 
theory with the assumption of large displacements (geometrical non-linearity) for the finite element simula- 
tion. Yet, this phenomenon does not have a strong effect on the final prediction of the flow-induced obstruction. 

The in-vitro setup allows to measure the pressure and the deformation during a fluid-structure interaction, in 
a configuration with strong asymmetry. Because of this asymmetry, the setup remains much more appropriate 
than other previous physical models, such as the collapsible tube (see e.g. [51129]). to study the human pharyn- 
geal airway. Furthermore, the measurement of the pressure at the level of the constriction, which plays a critical 
role in the flow-induced collapse, coupled to the measurement of the bidimensional geometry of the duct at the 
constriction, permitted to refine the results and analysis carried out in [6]. Nevertheless, as a simplification of 
the complex physiological reality, this setup has some limitations. For instance, it neglects the tridimensional 
effects that should be involved in a true human airway, where lateral walls seem to play a role during the 
collapse [3j. In the protocol, it would also be interesting to measure the lateral deformation of the tongue 
replica, so as to compare it to the numerical prediction. More pressure sensors would also help to capture with 
more accuracy the pressure distribution. 

Finally, concerning the clinical implications of this work, the numerical model proposed, as other models based 
on continuum mechanics (see e.g. |22j). is able to take into account easily and with relative accuracy patient- 
specific geometrical and mechanical properties. In particular, the choice of the Hookc law to model the soft 
tissue is reasonable for patient-specific modelling, since it appears difficult to obtain more information than the 
Young modulus from in- vivo measurements. A second advantage of the proposed model is that the simulation 
time is reduced due to the simplificatory assumptions. It is of the order of 20 min on a standard computei@, 
which may already be satisfying for a physiologist or a clinician. With a simple optimization of the code, still in 
a preliminary version, this time would be expected to be the order of one minute. At last, the in-vitro validation 
described in this paper should normally be an argument so that a clinician or a physiologist may trust the 
predictions from the numerical model, though in-vivo validation in this perspective is a must and remains our 
ultimate goal. 



5 Conclusion 

Extensive experimental validation of a numerical model that predicts the flow-induced collapse of the pharyn- 
geal airway in conditions of strong obstruction has been carried out. An in-vitro setup which reproduces the 
asymmetries and the particularities of the airway at the base of the tongue in pathological conditions has been 
used for this purpose. The prediction accuracy for the obstruction is of 4.2 % (mean quadratic error between 
prediction and measurements concerning the variation of the constriction area), and has been assessed for a 
wide range of parameter values. Then, it results from these comparisons that the numerical model may be 
considered as satisfying in a first approximation to predict the flow-induced deformation. Since it is based 
on simplified assumptions, a low computational cost is associated to each numerical simulation, which is an 
advantage for clinical applications. 
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Table 1 Comparison between the predictions and the measurements of the pressure at the constriction P c (a) and the 
constriction area A c (b), for different values of P ex t- 
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x=0 
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Fig. 1 Sagittal view of the upper airway overlaid with an outline of the in- vitro setup (a). Outline of the in- vitro setup, 
from frontal view (b). Experimental measurement of the wall deformation with a digital camera (c). The acquisitions 
of the camera are treated automatically with a software. Pressure sensors allow to obtain at each time step information 
about the airflow, which complements the information that comes from the camera. 
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(c) 

Fig. 2 Photographs of the in-vitro setup, (a) Deformable upper cylinder: view of the deformable portion, (b) The 
mounting of the 'in-vitro' tongue replica, (c) Frontal view of the replica at downstream end. 
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3 Finite element model of the latex wall of the in-vitro setup, with boundary conditions. The nodes indicated by 
arrows are immobilized, as they are supposed to be fixed to the rigid pipe. 
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